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INTRODUCTION 



At large distance scales one expects effectively: 



Confinement is wired in the structure of the ground state 
of the Yang-Mills (YM) theory. In the Schrodinger pic- 
ture in a physical gauge, the information is carried by the 
vacuum wave functional. One can approach the problem 
from various angles: 

1 . Try to develop a systematic strong-coupling expan- 
sion of the vacuum wave functional (VWF) [1]. 

2. Express the theory in cleverly selected variables and 
look for the VWF in terms of these variables [2]. 

3. Make a variational Ansatz for the wave functional 
in a certain gauge and find its parameters by minimizing 
(e.g.) the expectation value of the YM Hamiltonian [3]. 

4. Guess an approximate form of the wave functional 
and test its consequences [4, 5]. 

We follow the last avenue: we have proposed a guess 
of a simple vacuum wave functional in the temporal 
gauge [5], and will present here a few pieces of evidence 
in favor of the claim that this simple form is not far 
from the true vacuum wave functional of the Yang-Mills 
theory. 

In the hamiltonian formulation in D = 2 + 1 dimen- 
sions and temporal gauge, the vacuum wave functional 
satisfies the Schrodinger equation: 
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This form has the property of dimensional reduc- 
tion [6]: The computation of a spacelike loop in 
2 + 1 dimensions reduces to the calculation of a Wilson 
loop in Yang-Mills theory in 2 euclidean dimensions, 
which is known to exhibit the area law. However, the 
true VWF cannot be that simple, as it leads to incorrect 
physical predictions at short and intermediate distances. 



A GUESS AT AN APPROXIMATE VWF 

The starting point for our guesswork [for SU(2) gauge 
theory] has been the QED vacuum wave functional [7], 
whose simplest generalization reads 
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where V ab is an adjoint operator which fulfills, for van- 
ishing coupling, 
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where < 3f c \A\g\ is the covariant derivative in the adjoint 
representation, < 2) 2 = ^ ■ ^ the adjoint covariant lapla- 
cian, An is the lowest eigenvalue of (— & 2 ), and m is a 
constant (mass) parameter proportional to g 2 ~ 1 /j3 . The 
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FIGURE 1. Cumulative data for C0mc vs - P m physical 
units, on lattices of extensions L = 16,24,32,40,48, and eu- 
clidean lattice couplings /?£ = 6,9,12. The curve represents 
(Oqq (p) using the parameters of m and g 2 quoted in the text. 



expression for the VWF is written in the continuum no- 
tation, but assumed to be properly defined on a lattice. 

In Ref. [5] we have provided a series of analytical 
arguments in favor of the proposed form: 

1. By construction, ^oiA] becomes the VWF of elec- 
trodynamics in the free-field limit (g — > 0). 

2. It is a good approximation to the true vacuum also 
for strong fields constant in space, varying only in time. 

3. The part of the VWF that depends on slowly varying 
fields B s iow takes on the dimensional -reduction form. The 
fundamental string tension, at a given /3, is then easily 
computed as Op = 3m/4j3. 

4. If one takes the mass m in the VWF as a free 
variational parameter and computes (approximately) the 
expectation value of the YM hamiltonian, one finds that 
a non-zero (finite) value of m is energetically preferred. 



NUMERICAL EVIDENCE 

To estimate how good or bad the proposed approximate 
vacuum state is, we have also performed some numerical 
tests. We compare a set of quantities computed in two 
ensembles of lattice configurations: 

1 . "Recursion " lattices — independent 2-d lattice con- 
figurations generated with the probability distribution 
given by the proposed VWF, for a choice of its param- 
eters m and /3 = 4/g 2 . The recursion method was pro- 
posed and described in detail in Ref. [5]. 

2. Monte Carlo lattices — 2-d slices of configurations 
generated by MC simulations of the 3-d euclidean SU(2) 
lattice gauge theory with the standard Wilson action with 
coupling Pe \ from each configuration, only one (random) 
slice at fixed euclidean time was taken. 



To compare results for these two ensembles, one needs 
a way of fixing parameters (»J,/3) of the VWF that cor- 
respond to the Wilson-action coupling /?£ of the MC en- 
semble at a fixed lattice size L. The simplest way, used in 
Refs. [5, 8], is to choose j3 = Pe, and fix m at given j3f 
and L to get the correct value of the fundamental string 
tension Gf(Pe,L)- (This will be called variant A below.) 
With this choice we computed the equal-time connected 
B 2 -B 2 correlator and determined the value of the mass 
gap from a best fit to its exponential fall-off at large dis- 
tances. The result for recursion lattices was compared to 
the values of the + glueball mass computed in simula- 
tions of the 3-d YM theory by Meyer and Teper [9]. The 
deviations were at the level of at most 6%. 

Another possibility (variant E) is to find (m,j5) cor- 
responding to a given (L,Pe) by using numerically de- 
termined values of the true VWF for some trial gauge- 
field configurations (non-abelian constant fields, abelian 
or non-abelian plane waves). The square of the VWF can 
be computed numerically in simulations of the 3-d YM 
theory by the method proposed long ago by Greensite 
and Iwasaki [10]. Details of the method and results for 
various types of configurations will be presented else- 
where [11]. We will only use results for a set of abelian 
plane waves with the longest wavelength at a given lat- 
tice size L, A = L, and varying amplitudes: 
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The measured squared VWF at a given /3s can be 
parametrized as I^mc^'H 2 « exp(-ft M c with: 

Rmc[U U) ] = 2 (a + Y.i)(OMc(p) + const (8) 

with p 2 = 2 (l — cos ^5-), while for our proposed wave 
functional 

Rgo[U u) ] = 2(o + yj)taao(p) + const, (9) 
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Values of (m,g) can be determined by the best fit of the 
data for (Omc(p) at a given (L,j3g) by the function (10). 
The quality of the fit is illustrated in Fig. 1, the parame- 
ters of our wave functional in physical units come out to 



TABLE 1. Parameters and lattice sizes. 
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FIGURE 2. The Coulomb-gauge ghost propagator. 
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FIGURE 3. The color-Coulomb potential for k: = 10. 



be gp hys = 1.465,m p hy S = 0.771. In lattice units we then 

use: 

fl(j8E,L) = v/of(/3 E ,L)/0.44GeV, (11) 

m(PB,L) = m ph y S a(/3 E ,L), g 2 (PE,L) = gp hys a(/3 E ,L). 

The sets of parameters used in further tests are summa- 
rized in Table 1 . 

For both recursion and Monte Carlo lattices, we con- 
centrated on a subset of quantities in the Coulomb 
gauge [8]. As argued by Gribov [12] and Zwanziger [13], 
the low-lying spectrum of the Faddeev-Popov opera- 
tor in Coulomb gauge probes properties of non-abelian 
gauge fields crucial for the confinement mechanism. 
The ghost propagator in Coulomb gauge and the color- 
Coulomb potential are directly related to the inverse of 
the Faddeev-Popov operator, and play a role in various 
confinement scenarios. In particular, the color-Coulomb 
potential represents an upper bound on the physical po- 
tential between a static quark and antiquark [14]. 

Fig. 2 displays results for the ghost propagator in 
Coulomb gauge at j3 = 9 on 32 2 lattice. It was computed 
(in coordinate space) in the usual way from the inverse of 
the Faddeev-Popov operator (in the subspace orthogonal 
to trivial constant zero modes due to lattice periodicity). 
The agreement of the ghost propagator computed in both 
sets of lattices is almost perfect, for both variants of the 
choice of parameters of the approximate VWF 

The situation for the color-Coulomb potential is more 
complicated. There exist rare "exceptional" configura- 
tions with a very low (though still positive) lowest non- 
trivial eigenvalue of the Faddeev-Popov operator. These 
configurations were extremely difficult to gauge-fix to 
the Coulomb gauge. If one evaluates the potential in each 
single configuration, the exceptional ones possess a very 
high absolute value of the potential at the origin, |V(0)|. 
One can then classify configurations by their values of 
| V (0) | , and evaluate average potentials from sets of con- 
figurations satisfying a number of cuts |V(0)| < K. Fig. 3 



shows results for K = 10 (satisfied by about 80% lattices) 
at j3 = 9 (32 2 lattice). The potentials agree quite well. 
However, the agreement between recursion and Monte 
Carlo lattices deteriorates with increasing K. Still, the po- 
tentials stay almost identical for recursion lattices with 
both variant choices of VWF parameters. 

In conclusion, the proposed vacuum wave functional 
for the temporal-gauge SU(2) YM theory in 2+1 dimen- 
sions, Eqs. (4, 6), is a fairly good approximation to the 
true ground state of the theory. We have accumulated an- 
alytical arguments and numerical results in its favor. The 
agreement of Coulomb-gauge quantities in recursion and 
Monte Carlo lattice ensembles is quite satisfactory for a 
bulk of the probability distribution, but there seems to be 
some disagreement in its tail. 
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